Submitted to ApJ 

Preprint typeset using I^T^^ style cmulateapj v. 10/09/06 



O 

o 

(N 

oo 
(N 

> 

O 

m 

00 

o 
o 

Oh 

6 ■ 



CASCADING OF FAST-MODE BALANCED AND IMBALANCED TURBULENCE 
T. K. Suzuki^, A. Lazarian^, A. Beresnyak^ 

Submitted to ApJ 

ABSTRACT 

We study the cascading of fast MHD modes in magnetically dominated plasma by performing 
one-dimensional (ID) dynamical simulations. We find that the cascading becomes more efficient as 
an angle between wave vector and underlying magnetic field increases and fast mode becomes more 
compressive. We also study imbalanced turbulence, in which wave packets propagating in one-direction 
have more power than those in the opposite direction. Unlike imbalanced Alfvenic turbulence, the 
imbalanced fast mode turbulence shows faster cascading as the degree of imbalance increases. We 
find that the spectral index of the velocity and magnetic field, which are carried by the fast-mode 
turbulence, quickly reaches stationary value of —2. Thus we conclude that the dissipation of fast 
mode, at least in ID case, happens not due to weak or strong turbulent cascading, but mostly due to 
nonlinear steepening. The density fluctuations, which are carried by slow-mode perturbation in the 
larger scale and entropy-mode perturbation in the smaller scale, depend on initial driving spectrum 
and a ratio of specific heat. 

Subject headings: magnetohydrodynamics - plasma - turbulence 



L INTRODUCTION 

Magnetohydrodynamical (MHD) turbulence is ubiq- 
uitous in astrophysics. All interstellar medium (ISM) 
phases are well magnetized with Larmor radius of ther- 
mal proton much smaller than outer scale. Molecular 
clouds are no exception, with compressible, often high- 
Mach number, magnetic turbulence determining most 
of their properties (see Elmegreen & Falgarone 1996, 
Stutzki 2001). Star formation (see McKee & Tan 2002; 
Elmegreen 2002; Pudritz 2002; Ballesteros-Paredes et 
al.2006), cloud chemistry (see Falgarone 1999 and ref- 
erences therein) , shattering and coagulation of dust (see 
Lazarian & Yan 2002 and references therein) are exam- 
ples of processes for which knowledge of turbulence is 
absolutely essential. It is also believed that MHD waves 
and turbulence are important for the acceleration of the 
solar wind (e.g. Tu & Marsch 1995) as well as winds from 
stars which possess surface convective layers (e.g. Tsuji 



As MHD turbulence is an extremely complex process, 
all theoretical constructions should be tested thoroughly. 
Until very recently matching of observations with theo- 
retical constructions used to be the only way of test- 
ing. Indeed, it is very dangerous to do MHD turbulence 
testing without fully 3D MHD simulations with a dis- 
tinct inertial range. Theoretical advances related to the 
anisotropy of MHD turbulence and its scalings (see She- 
balin et al. 1983, Higdon 1984, Montgomery, Brawn & 
Matthaeus 1987, Shebalin & Montgomery 1988, Zank & 
Matthaeus 1992) were mostly done in relation with the 
observations of fiuctuations in solar wind. Computers al- 
lowed an important alternative way of dealing with the 
problem. While still presenting a limited inertial range, 
they allow to control the input parameters making it eas- 
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ier to test theoretical ideas. 

It is well known that linear MHD perturbations can 
be decomposed into Alfvenic, slow, fast, and entropy 
modes with well-defined dispersion relations (see Alfvcn 
& Falthammar 1963). The separation into Alfven and 
pseudo-Alfven modes, which are the incompressible limit 
of Alfven and slow MHD modes, is an essential element 
of the Goldreich-Sridhar (1995, henceforth GS95) model 
of turbulence. There the arguments were provided in fa- 
vor of Alfven modes developing a cascade on their own, 
while the pseudo-Alfven modes being passively advected 
by the cascade. The drain of Alfvenic mode energy to 
pseudo-Alfven modes was shown to be marginal. 

The separation of MHD perturbations into modes in 
compressible media was discussed further in Lithwick & 
Goldreich (2001) and Cho & Lazarian (2002, 2003 hence- 
forth CL02, CL03, respectively). Even though MHD tur- 
bulence is a highly non-linear phenomenon, the modes 
does not constitute an entangled inseparable mess. 

The actual decomposition of MHD turbulence into 
each mode was a challenge that was addressed in CL02, 
CL03. They studied a particular realization of turbu- 
lence with mean magnetic field comparable to the fluctu- 
ating magnetic field. This setting is, on one hand, rather 
typical of the most MHD flows in Galaxy and, on the 
other hand, allows them to use a statistical procedure of 
decomposition in the Fourier space, where the bases of 
the Alfven, slow, fast, and entropy perturbations were 
defined. The entirely different way to decompose modes 
in high or low-/? cases in real space was shown in CL03 
to correspond well to this procedure'^. 

In particular, calculations in CL03 demonstrated that 
fast modes are marginally affected by Alfven modes and 
follow acoustic cascade in both high and low f3 medium. 
In Yan & Lazarian (2002; 2004; henceforth YL04) the 
fast modes were identified as the major agent in scatter- 
ing of cosmic rays. Similar results in relation to stochas- 

^ In CL03, an isothermal equation of state is used so that 
entropy-mode fluctuation is prohibited. 
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tic acceleration of cosmic rays were reached in Cho & 
Lazarian (2006). Interestingly enough, the dominance of 
fast modes for both acceleration and scattering is present 
in spite of the dissipative character of fast modes. The 
reason for that is isotropy** of this mode reported in 
CL03. The Alfven modes, that are still the default for 
scattering and acceleration for many researchers, are in- 
efficient due to the elongation of eddies along the mag- 
netic field (Chandran 2000, Yan & Lazarian 2002). 

Apart from cosmic rays, acceleration of charged inter- 
stellar dust is also dominated by fast modes (Lazarian 
& Yan 2002, Yan & Lazarian 2003). In stellar magneto- 
spheres with spiral magnetic fields, fast-modes possibly 
play a role in energetics and dynamics of winds from 
usual stars ("Suzuk i et al.ll2006l ) and pulsars (|Lvubarskv I 
[2003). These issues provide an additional motivation to 
studies of properties of fast mode cascade. 

The most efficient interaction of fast modes is expected 
to happen for waves moving in the same direction (see 
CL02). This indicates that, unlike studies of Alfven 
modes for which the 3D character of interactions is es- 
sential, ID simulations could still give, at least in some 
respects, correct physical properties of the fast-mode tur- 
bulence. ID simulations have a big advantage that we 
can use a number of grid points for a broader inertial 
range than 2D or 3D simulations. Based on these issues, 
we carry out ID decay simulations of fast-mode turbu- 
lences. We should cautiously analyze simulation data, 
because the ID geometry systematically enhances shocks 
without dilution to the tangential directions. 

In what follows we discuss our code in §2, calculate 
power spectra in §3, and study cascading of fast modes 
in §4. We discuss our findings in §5. 



2. THE CODE AND THE SIMULATIONS SETUP 

We initially give perturbations of fast mode by fluc- 
tuations of longitudinal velocity, Vx , transverse velocity, 



dp 



dv,. 



transverse magnetic field. By, density, p, and pres- 
sure p (or specific energy, e = ^^^j^^^ , where 7 is a ra- 
tio of specific heat) in magnetically dominated plasma, 
namely low-/? conditions. We dynamically treat time 
evolution of each variables by solving ideal MHD equa- 
tions without continuous driving (decay simulation) . We 
do not consider the third (z) component of magnetic 
field and velocity, which indicates that Alfven mode is 
switched off in our simulation, whereas fast, slow, and 
entropy modes are automatically treated. We assume 
one-dimensional approximation, namely all the physical 
variables depend only on x; a wave number vector has 
only the x-component, namely, k = k^- The number of 
grid points of the x component is 65536(= 2^^), and the 
periodic boundary condition is imposed in x-dircction. 

* Recent calculations of fast modes interacting with Alfven 
modes by Chandran (2005) indicate some anisotropy of fast modes. 
In particular, he shows that contours of isocorrelation for fast 
modes get oblate with shorter direction along the magnetic field. 
However, the calculations are performed assuming that Alfven 
modes are in the weak turbulence regime. This regime has a lim- 
ited inertial range. Moreover, anisotropic damping discussed in 
YL04 is probably a more strong source of anisotropy. Anyhow, 
the above conclusions about the dominance of fast modes are not 
altered by the anisotropies that corresponds to squashing contours 
of isocorrelations in the direction of magnetic field. 
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where B"^ — B^.^ By and v"^ — 
Lagrangian and Eulerian time derivatives. 

The minimum wave number corresponds to the box 
size. We set-up p = 1, p — 1, Vx — 0, Vy = a,s back- 
ground conditions. We test two cases of the background 
magnetic fields, (Bx,By) = (9,3) (quasi-parallel) and 
{Bx,By) = (3,9) (quasi-perpendicular), and three cases 
of a ratio of specific heat, 7 = 5/3 (adiabatic), 1.1 (nearly 
isothermal), and 1 (isothermal). 

For the simulations we employed 2nd order nonlinear 
MHD Godunov-MOCCT scheme developed by Sano & 
Inutsuka (2007, in preparation) (see also Suzuki & Inut- 
suka 2005; 2006). In this scheme, each cell boundary is 
treated as discontinuity, for the time evolution we solve 
nonlinear Riemann shock tube problems with magnetic 
pressure by the Rankine-Hugoniot relations. Therefore, 
entropy generation, namely heating, is automatically cal- 
culated from the shock jump condition. A great advan- 
tage of our code is that no artificial viscosity is required 
even for strong MHD shocks; numerical diffusion is sup- 
pressed to the minimum level for adopted numerical res- 
olution. 

3. TIME EVOLUTION AND POWER SPECTRA 

We firstly study time-evolution of the simulated tur- 
bulences and their power spectra. In this section, we 
present the results of quasi-parallel {{Bx,By) — (9,3)) 
and 7 = 1.1 case. The initial amplitude is set to 
be dBy — 4.4, corresponding to super-sonic and sub- 
Alfvenic condition. Fast modes traveling in both di- 
rections have equal energies, we call this balanced tur- 
bulence. We test two types of initial spectral slopes; 
one is oc k~^/'^ with random (5-correlated phases for all 
wavenumbers, and the other is oc fc*^, namely white noise. 
Note that /c""^/^ is a prediction of weak acoustic turbu- 
lence, which is supposed to have similar propertiles to 
weak fast-mode turbulence (CL02; CL03), whereas the 
validity of weak turbulence appro ximation is still un- 
der debate for acoustic turbulence (iFalkovich &: Mever I 



119961: IZakharov. L'vov. fc Falkovich ||1992D . Figures [Tjfc 
[2] present time-evolution of the cases with initial spectra 
oc fc~^/^ and k^ , respectively; the six panels present Vx, 
Vy, By, p, e, and p at t = 0.5 (dashed) and t = 3.0(solid). 

We also show the power spectra of v, B, and p at 
t = 3.0 in Figures [3] and S) Although our initial con- 
ditions contain only fast mode, slow mode is generated 
by nonlinear process of mode interaction. In order to un- 
derstand the interplay of modes, we decompose all val- 
ues into fast and slow perturbations (see CL02). This 
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Fig. 1. — Evolution of Dj; (top-left), Vy (middle- left), By (bottom- 
left), p (top-right), e (middle-right), and p (bottom-right). Dashed 
and solid lines are the results at t = 0.5 and t = 3, respectively. 
The initial spectral slope is oc k"^!"^ . 



0.05 









^ v\ 


r - 


0.05 


>\ /vi / A / 

-1; ■'i I, 

: 1/ 


\ 1 /" 



0.2 






0.1 

>°~o 






-0.1 






-0.2 









0.05 0.1 0.15 0.2 0.1 



Fig. 2. — The same as Figure [2] but for white noise as the initial 
perturbation. Note that the horizontal axis is enlarged compared 
with Figure|2] the horizontal axis of the right panels is from a:: = 0.1 
to 0.2, while that of the left panels is from x = to 0.2. 

method is statistically correct even for nonlinear waves 
as discussed in CL03. In Figures[3]and[4l the decomposed 
spectra of fast (left) and slow (middle) modes as well as 
the raw spectra (right). Note that the By and p spectra 
of the slow mode in the initial white noise case (Figure[4]) 



may contain errors in the small scale (large k) region. As 
we describe later, the density perturbation is mainly due 
to entropy-mode in this case. The slow mode spectrum 
is calculated from the substraction of the dominated en- 
tropy mode from the total power, which might leave an 
error in the slow mode. The By spectrum of the slow 
mode, which is connected to the slow mode density spec- 
trum through the frozen-in condition, might also contain 
non- negligible errors. 

The velocities and magnetic fields exhibit many discon- 
tinuities, i.e. fast and slow MHD shocks, in both cases 
(the left panels of Figures [T] & [2|) . The amplitudes of the 
white noise case is much smaller than those of the fc~^/^ 
spectrum case, because the white noise perturbations ini- 
tially contain more energy in smaller scales which suffer 
faster damping. 

Figures [3] & [4] indicates that the velocity and mag- 
netic field perturbations are carried by fast mode. This 
is first because the slow mode is downconverted from the 
fast mode, and second because the fast wave becomes 
magnetic mode in low-/3 medium. The slow wave es- 
sentially corresponds to hydrodynamic (acoustic) mode, 
and this is the main reason why the slow-mode pertur- 
bation dominates the fast-mode in the density spectra, 
in accordance with the earlier claims in CL03, Passot & 
Vazquez- Semadeni (2003), and Beresnyak et al. (2005), 
whereas the entropy-mode perturbation is also important 
in the initial white noise case as discussed later. 

As expected from many discontinuities in Figures [T] & 
[21 V and By exhibit shock dominated spectra, cx in 
both fc~^/^ (Figure 131) and white noise (Figured]) cases. 
This indicates that the weak cascade solution of acous- 
tic turbulence (cx fc~^/^) does not hold for any significant 
time even though it is provided initially. The randomness 
of phases typically assumed in acoustic turbulence or in 
other theories of weak turbulence breaks down quickly 
by nonlinear interaction, leading to Burgers-like shock- 
dominated random state. This is different from the re- 
sults of the 3D simulation by CL02 in which spectrum 
fc~^/^ is suggested. The difference may be due to the 
difference of ID vs. 3D cases; the effect of shocks tend 
to become systematically predominant in ID simulations 
because shocks do not dilute geometrically. 

Since our simulations do not assume isothermal gas, 
the different dissipation characters influence the evolu- 
tion of thermal properties. The initial specific energy. 



eo(= 



), is 10 in both cases. In the k '^/^ case, the 



(7-l)pQ ' 

heating is localized at first in places where strong shock 
heating takes place; at t = 0.5, e partially becomes > 12, 
while some places stay e w eo(= 10). Most places expe- 
rience the heating eventually and at t = 3.0 the gas is 
averagely heated up to e « 12. On the other hand, in the 
white noise case (Figure [5]), e is quickly heated to e « 17 
becauase of the rapid dissipation, and after that e stays 
almost constant with time. 

The densities show different features in Figures[l]and[2l 
the p structure of the initial k~^/'^ spectrum case (Figure 
[T]) is dominated by larger scale with a couple of shocks, 
which is qualitatively simillar to v and By, while many 
small-scale structures are predominant in the white noise 
case (Figure[2]). Figure [3] (fc~'^/^ case) shows the density 
perturbation is mainly carried by the slow-mode and the 
spectral slope is the shock-dominated one, k~'^ , similarly 
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Fig. 3. — Power spectra of v (top), B (middle), and p (bottom) 
decomposed into fast (left) and slow (middle) modes at t = 3.0 for 
the intial oc k~^^^ perturbation. Composed values are shown on 
the right for comparison. 
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Fig. 4. — Same as Figure[3]but for the initial white noise pertur- 
bation. Nte that the composed density spetrum is mainly owing to 
entropy mode (note show); this is the reason why the total power 
of the density turbulence is larger than the combination of the fast- 
and slow-mode perturbations. The slow mode By and p spectra 
might contain errors in the high k region (see text). 
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Fig. 5. — Density spectra of the initial white noise perturbations 
for 7 = 1.1 (dotted) and 1 (dashed) at t = 3. 

to the By and v spectra. 

On the other hand, in the initial white noise case, the 
total power of the density fluctuation (the right bottom 
panel of Figure 2]) is much larger than the sum of the 
fast- and slow-mode perturbations. In fact, the density 
perturbation in this case is mostly owing to an entropy 
mode, which is a zero-frequency mode with unperturbed 
p and with perturbed p and e (see e.g. Lithwick & Gol- 
dreich 2001). By a close look at Figure one can see 
that the phases of e and p, both of which show small- 
scale structure, are anti-correlated so that p(oc pe) shows 
a smooth feature. The initial perturbations (white noise) 
contain sizable energy in the small scale. These small- 
scale turbulences quickly dissipate, rather than converted 
to slow-mode perturbation, to entropy-mode fluctuation. 
The density fluctuation of the entropy mode drives the 
fluctuation of the energy deposition by shocks. The heat- 
ing rate per unit mass is larger in hotter regions because 
of the lower density, which self-sustains the entropy-mode 
perturbation. As a result, the density perturbation is 
preserved with keeping the initial flat spectrum until the 
end of simulation (t = 3), when the turbulent energy of 
the fast mode decays ~ 1/1000 of the initial value. 

In conclusion, the density perturbation is controlled 
by initial driving turbulence (see Beresnyak 2005 for 3D 
case). This is a kind of bistability. When the energy is 
given in a larger scale, the density perturbation is carried 
by the slow-mode downconverted from the fast-mode. 
The generated slow-mode shows the shock-dominated 
spectrum. On the other hand, when the energy is in- 
put in a smaller scale, the density perturbation is self- 
sustained by the entropy-mode fluctuation, which is a re- 
sult of the dissipation of the fast-mode turbulence. The 
density perturbation triggers the fluctuation of the heat- 
ing, which sustains the entropy-mode fluctuation in it- 
self. In this case, an initial spectrum is preserved in the 
density perturbation. 

In isothermal (7 = 1) gas, the entropy mode does 
not exist because p must follow unperturbed p; the p 
spectrum would be different especially in a small scale. 
We carry out the simulation of isothermal gas with the 
initial white noise fluctuation. Figure [5] compares the 
density power spectra adopting the initial white noise 
perturbations for 7 — 1 (isothermal) and 1.1 at t = 3. 
This figure shows that the density spectrum, which is 
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Fig. 6. — Energy of the fast modes of various initial amplitudes. 
The background field strength is {Bx,By) = (9,3). Dot-dashed, 
solid, and dashed lines are results of dBy = 0.5 (sub-sonic), 4.4 
(super-sonic/sub-Alfvenic ), and 13.2 (super- Alfvenic ). The initial 
energy density is normalized at unity in each case. 

due to the slow-mode in the isothermal gas, becomes 
steep nearly cx fc^^. Therefore, the imprint of the ini- 
tial spectrum is unique in the density perturbations of 
non-isothermal gas, in which entropy-mode perturbation 
does exist. We should note, however, that the fast-mode 
turbulence is not influenced whether we adopt isothermal 
or nonisothermal conditions. 

The above result shows that density perturbation in a 
small scale is influenced by an equation of state of gas. 
In general, a ratio of specific heat, 7, is controlled by 
physical processes such as radiative cooling and thermal 
conduction. This shows that we might get some infor- 
mation on gas properties from observations of density 
turbulence. Density perturbation is also important in 
parametric instability because it can work as a mirror 
against Alfven and magnetosonic waves. Thus, the re- 
lation between 7 and density turbulence is interesting. 
However, unlike fast-mode turbulence, 3D treatment is 
essential to entropy as well as slow modes, since they are 
passively advected (jLithwick fc Goldreich |[2"00 1^ and be- 
have very differently in ID and 3D; we pursue this issue 
in our future project. In this paper we concentrate on 
the evolution of the fast mode from now. 

4. CASCADING OF FAST MODE 

We analyze time-evolution of integrated energy density 
(energy column density), 

J dx{pSv^ + 5B'^)/2, of the decaying fast modes. 

4.1. Dependence on Initial Amplitude 

Figure [6] shows dissipation of fast mode energy as a 
function of time for different initial amplitudes with the 
same background {Bx,By) = (9,3). The figure indicates 
that larger amplitude gives faster damping. This is rea- 
sonable since the dissipation is owing to nonlinear pro- 
cesses such as shocks and wave- wave interactions. Damp- 
ing times, which are defined as slopes of the respective 
lines, become longer at later times as the amplitudes de- 
crease, which is also consistent with the nonlinear dissi- 
pation through steepening. 



Fig. 7. — Energy (per grid) of the fast modes of the quasi-parallel, 
{Bx,By) = (9,3), (solid) and quasi-perpendicular, {B^^By) = 
(3,9), (dashed) cases. 

7 = 1.1 are adopted. The figure shows that the dissi- 
pation is more effective in the quasi-perpendicular case. 
This is because the fast mode shows more compressive 
character as an angle between the propagation and the 
magnetic field increases and the shock dissipation is en- 
hanced (Suzuki et al. 2006). 

4.3. Dependence on Plasma-(3 

Since our simulations do not assume isothermal gas, 
the gas is influenced by the heating, and the plasma- 
P values change with time by the energy transfer from 
magnetic fields to gas. This effect is more prominent 
in the adiabatic cases (7 = 5/3). We have found from 
the simulations that the dependence of the decay of the 
fast-mode turbulence on plasma-/? is very weak provided 
the plasma is kept magnetically dominated {j3 < 0.5) 
during simulations. However, if (3 approaches unity, the 
dissipation becomes saturated. 

Figure [8] compares cases with the same {Bx,By) — 
(9,3)(left) or {B^,By) = (3,9)(right) and dBy = 13.2 
but different 7 = 1.1 (dotted) and 5/3 (solid). In the 
case with 7 = 5/3 the plasma is heated up so that the gas 
pressure eventually becomes comparable with the mag- 
netic pressure (/? ~ 1). In contrast, the plasma (3 stays 
well below unity in the cases with 7 = 1.1 because the 
heating is less efficient due to the smaller 7. All the 
three panels of Figure [8] show that once the plasma-/? 
becomes unity the dissipation is suppressed in the case 
with 7 = 5/3, while the energy density monotonically de- 
creases in the case with 7 = 1.1. In the middle and right 
panels, the fast-mode energies of the 7 = 1.1 cases are far 
below those of the 7 = 5/3 cases. Even in the left panel, 
the energy of the 7 = 1.1 case would be smaller than 
that of the 7 = 5/3 case, if we proceeded the simulation 
longer time. 

At later epochs in the case with 7 = 5/3, slow mode 
also has a sizable amount of energy because the phase 
speed of the slow mode becomes comparable with that 
of the fast mode. This tendency is more extreme in the 
quasi-perpendicular case (the right panel), because the 
plasma is heated up more due to larger compressibility. 

4.4. Balanced vs. Imbalanced 



4.2. quasi- Parallel vs. quasi-Perpendicular 

Figure [7] compares the energy of the fast modes of 
quasi-parallel (B^ — 9, By = 3) and quasi-perpendicular 

{Bx — i^By — 9) cases. The same dBy — 4.4 and total fast mode energy in balanced case with dBy = 4.4 



Figure [9] compares the evolution of fast modes in bal- 
anced and fuUy-imbalanced cases. The background field 
strength is {Bx,By) — (9,3). Dot-dashed line indicates 
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Fig. 8. — left : Energy (per grid) of the fast modes of different 7 = 
1.1 (dotted) and 5/3 (solid). The same background field strength, 
{Bx, By) = (9, 3), amplitude, dBy = 13.2, and the initial spectrum, 
oc fc"'^/^, are adopted, middle : the same as the left panel but for 
{Bx,By) = (3,9). right : the same as the left panel but for the 
initial white noise spectrum. 
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Fig. 9. — Comparison of balanced and imbalanced fast-mode 
turbulence. Dot-dashed line indicates total fast mode energy in 
balanced case with dBy = 4.4, solid line indicates fast mode energy 
propagating to one direction in balanced case with dBy = 6.2, 
and dashed line denotes fast mode energy in imbalanced case with 
dBy = 4.4. 
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Fig. 10. — Power spectra of the balanced (solid) and balanced 
(dotted) turbulence at t = 3. 

solid line indicates fast mode energy propagating to one 
direction in balanced case with dBy ~ 6.2, and dashed 
line denotes fast mode energy in imbalanced case with 
dBy = 4.4. Note that the energy is proportional to dBy 
so that the initial values of these three lines are the same. 
The figure shows that the imbalanced cascade is more 
dissipative, which is opposed to the tendency obtained 
in the Alfvenic turbulences. 

One reason is different nature of the dissipation of fast 
and Alfven modes. The fast mode is isotropic. The dissi- 
pation occurs among fast waves which meet the following 
resonance conditions (see CL02): 

Wl+CJ2=C^3, (6) 
fcl+fc2=fc3- (7) 

Since, uj k for the fast modes, the resonance con- 
ditions can be satisfied only when all three k vectors 



are collinear. When parent '1' and '2' waves propagate 
in the counter directions, (different signs of ki and k2) 
the daughter wave should have smaller wavenumber, k^, 
(larger wavelength). 

In contrast, if parent waves propagate in the same 
direction, wavenumber of the daughter wave becomes 
larger than the parents', which is more appropriate for 
the dissipation. In the imbalanced case, we only have 
fast modes traveling in the same direction so that the 
resonances selectively generate waves with shorter wave- 
lengths which suffer faster dissipation. Dot-dashed (en- 
ergy of one direction in balanced case with dBy = 6.2) 
and dashed (imbalanced case with dBy — 4.4) lines show 
very similar damping at the beginning [t < 0.5). This in- 
dicates that the dissipation of ID fast mode turbulence 
is controlled by wave packets traveling in the same di- 
rection and almost independent of those in the counter 
directions during this period, since the initial energy of 
one direction in the former case is identical to the initial 
total energy of the latter. 

In later time, however, the dissipation of the balanced 
cases is slower than the dissipation of the imbalanced 
case. To investigate this difference, we compare the 
power spectra at i = 3 (Figure [T(I| . The comparison 
shows that more power is kept at very small k « 10^^ in 
the balanced turbulence, although both show the shock- 
dominated spectra, cx in the higher k regime. This 
might indicate that a kind of inverse cascade takes place 
preferentially in the balanced turbulence. Since turbu- 
lences with larger scales (smaller k) suffer less damping, 
further dissipation seems suppressed in the balanced tur- 
bulence. However, we think that to discuss this issue the 
ID simulation is not sufficient, hence, we leave it to our 
future studies. 

We also suspect that bulk acceleration affects the faster 
dissipation of the imbalanced fast-mode turbulence. The 
imbalanced flux of fast mode transfers the momentum 
flux to the fluid. As a result, the system is accelerated 
into the direction corresponding to the initial wave mo- 
mentum flux. Generally, energy density of waves is lost 
in accelerating fluid by transfer of the momentum flux to 
the gas e ven though the waves themselves do not suffer 
damping (lJacqueslll977t l. However, after comparison of 
the bulk momentum flux of the gas with the dissipated 
energy of the fast-mode turbulence, we have found that 
the effect of the bulk acceleration is not so large (< 10%) 
in our case. 

5. SUMMARY AND DISCUSSIONS 

We have shown that the imbalanced fast modes dis- 
sipate in a different way than Alfven ones. For the 
Alfven turbulence cascading is due to wave packets mov- 
ing in the opposite direction, therefore, the imbalance 
makes the cascading less efficient. For fast modes waves 
moving in the same direction interact more efficiently. 
We confirmed that the oppositely moving fast modes are 
not essential for the cascading as discussed earlier (see 
CL02). 

However, the obtained power spectrum (fc^^) is dif- 
ferent from the weak acoustic turbulence (fc^"^/^) sug- 
gested in CL02, which may be owing to the different 
geometries (ID vs. 3D) as discussed in fJHl We should 
note that this problem of the scaling of acoustic turbu- 
lence is not settled down because the approximation of 
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weak cascade is controversial ( 


Falkovich & Mever II1996I: 


IZakharov. L'vov. & Falkovich 


119921). It is clear that fur- 



ther research by 3D simulations with a sufficient inertial 
range is necessary. 



We have also found that as the amplitude of com- 
pressions increase, the non-linear damping of turbulence 
speeds up. In our adiabatic calculations the restoring 
force was increasing due to medium heating and as the 
result the non-linear steepening and cascading slowed 
down. Similarly, for the modes that compress both mag- 
netic field and gas, the rate of dissipation was less com- 
pared to the modes that mostly compress the gas. 

The density perturbation is mainly carried by slow 
mode perturbation in the large scale and entropy mode 
in the small scale perturbation in non-isothermal gas. 
The slow-mode density spectrum is dominated by shocks, 
oc fc^^, while the entropy-mode, which is a consequence 
of the dissipation of the initial fast-mode turbulence, pre- 
serves the initial turbulent spectrum. As a result, the 
spectral index of the density perturbation is subject to 
the initial spectrum of given fast-mode turbulence. In 
isothermal gas, density spectrum is different in a small 
scale because entropy- mode fluctuation does not exist. 
This shows that to study density turbulence one has 
to carefully treat driving properties of turbulence and 
to take into account appropriate thermal processes (e.g. 
Koyama & lutsuka 2002; Larson 2005) rather than as- 
sume a polytropic equation of state. While the density 
turbulence is an interesting issue, obviously 3D treatment 
is required since the slow- and entropy-mode perturba- 
tion is different in ID and 3D circumstances; we will 
study this topic in a future paper. 

In this paper we have focused on the fast-mode in low- 
j3 medium. Various interstellar/intergalactic medium is 
dominated by magnetic pressure rather than gas pres- 
sure, i.e., low- [3 condition. Referring to the tabulation in 
YL04 (Table 1 in their paper), galactic halo, and warm 
and cold ISM are mildly low-/3 (/3 = 0.1 — 1), and dark 
clouds are more magnetically dominated, (i — 0.01 — 0.1. 
Our simulations are directly applicable to these media. 

Magnetospheres of general stars with surface convec- 
tion (like the sun) and compact objects (white dwarfs and 
neutron stars) are also in low-/3 condition, because the 
density rapidly decreases upwards. Moreover, the turbu- 
lence tends to be imbalanced there because it is mainly 
generated from central stars and dominated by the out- 
ward component. Our results show that the fast mode is 
likely to be the culprit in the heating and acceleration of 
winds because it dissipates effectively even without the 
inward component. The advantage of the fast-mode ver- 
sus Alfvenic turbulence in the stellar wind acceleration is 
that the latter requires production of counter-waves via 
reflection or by some other mechanism. 

Cascading of imbalanced fast modes is important for 
the problem of streaming instability evolution, which is 
the part of the picture of both galactic leaky box model 
and the acceleration of particles in shocks. This insta- 
bility arises from the flow of particles along the mag- 
netic field direction and reflects the particles back (Kul- 
srud & Pearce 1968). Our results indicate that the fast 
mode cascading limits the instability even in the absence 
of both damping and the ambient turbulence (YL02, 
Farmer & Goldreich 2004). 

In our simulation the dissipation is mainly due to MHD 



shocks besides sub-grid scale damping. However, if we 
discuss the dissipation in more detail, we need to take 
into account the microphysics. Galactic halo and stellar 
winds generally consist of coUisionless plasma. In such 
conditions, the fast-mode as well as slow-mode fluctu- 
ations possibly dissipate by transit-time damping (e.g. 
Barnes 1966; Suzuki et al.2006). In warm ionized ISM, 
the dissipation is mainly due to Coulomb collisions. If 
there is a considerable fraction of neutral atoms, e.g. in 
cold ISM and dark clouds, collisions between neutrals 
and ions dominate (see Table 1 in YL04) . 

Turbulence may be very different in ID and 3D. We 
find, however, that our attempts to get insight into the 
dynamics of fast modes using ID model is meaningful, as 
the possible transversal deviation 5k is limited. Assum- 
ing that the fast modes follow the acoustic cascade with 
the cascading time 

Tk = Lo/k\l = (fc/L)-i/2 X VphlV\ (8) 

YL04 used the uncertainty condition Scutcas ^ 1, where 
Slo ^ Vph5k{5k / k) to estimate 5k 

1 f V 

Eq. ^ provides a very rough estimate of how much 
the randomization of vectors is expected due to the fast 
mode cascading. This estimates suggest that deep down 
the cascade where kL ^ 1 the approximation of waves 
moving in the same direction could be accurate enough. 

In this paper we dealt with fast modes irrespectively 
of the evolution of the other modes. This possibil- 
ity corresponds to the earlier theoretical and numeri- 
cal studies (GS95, Lithwick & Goldreich 2001, CL02, 
CL03). Recently, however, Chandran (2005) discussed 
the interaction between fast and Alfven modes under 
the approximation that the latter constitute weak tur- 
bulence, i.e. turbulence that cascades slowly due to 
weak non-linearities. We expect appreciable mitigation 
of the Alfven-fast mode interaction compared to those 
in Chandran (2005), when Alfvenic turbulence is strong, 
i.e. when the turbulence evolves fast, over just one wave 
period. We feel that the resulting large disparity in the 
rate of evolution of the Alfvenic and fast modes should re- 
sult in much less cascading of fast modes by the Alfvenic 
modes compared to the case of the weak Alfvenic tur- 
bulence. Note, that strong Alfvenic turbulence is default 
for most astrophysical situations, while the inertial range 
of the weak Alfvenic turbulence is limited^ (Goldreich & 
Sridhar 1997, Lazarian & Vishniac 1999). In our next 
paper we shall provide calculations of the fast modes in- 
teraction with the Alfvenic modes. 
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